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Abstract 

In this paper, way-point tracking control of a container ship 
based on LOS method using adaptive stochastic sliding 
mode controlhas beeninvestigated.Effective control of ships 
in a designed trajectory is always an important task for ship 
maneuvering. The design is based on a low and high 
frequency model of the vessel motion adequate to ship 
steering. The low frequency model describes the vessel 
response to rudder control and slowly varying 
environmental forces. The high frequency model represents 
the wave induced oscillatory part of the yaw motion.For 
suppression of high frequency wave disturbances, the 
nonlinear recursive filter is considered. The uncertainty of 
model parameters in the presence of disturbances such as 
waves, wind and currents, dictates the application of 
techniques in which the nonlinear equations of ship's motion 
and the presence of unknown parametershave been taken 
into consideration. In this sense, the adaptive stochastic 
sliding modedesign is used for solving this problem. The 
simulation results show the suitability of this technique to 
compensate the time-varying disturbances effects in way- 
point tracking control. 
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Introduction 

The way-point tracking control problem is an issue of 
high interest in research areas of ship maneuverings. 
In seaway, the path or trajectory is constructed by 
using a set of way-points that can be generated 
according to sail plan and weather data (minimum 
resistance or energy approach). The way -point 
tracking control problem is to make the ship follow the 
path planned with the way-points by controlling the 
rudder. A widely used method for way-point tracking 
control is line of sight (LOS) guidance. In this 
methodology, a LOS vector is computed from the ship 
position to the next way point for heading control 
[Fossen (2000)]. The way-point tracking control has 


been improved over several years. Reference [Peterson 
& Lefeber (2001)] gives a yaw torque control law for 
ship way-point tracking control problem based on a 
full state feedback control approach. A new fuzzy 
autopilot for way-point tracking has beenproposed in 
[Cheng et al. (2006)]. A way-point guidance algorithm 
by LOS calculates a dynamic LOS vector norm in 
order to minimize the cross track error, is presented in 
[Moreira et al. (2007)]. In reference [Aguiar and 
Pascoal (2007)] integrator backstepping and Lyapunov 
based techniques are used for way-point tracking of 
under actuated autonomous underwater vehicles 
(AUVs) in the presence of constant unknown ocean 
currents and parametric modeling uncertainty. Model 
predictive control (MPC) for a way-point tracking 
under actuated surface vessels with input constraints 
wasproposed in [Ryeko, Jing (2010)]. In order for the 
control action to render good helmsman behavior, a 
MPC scheme with line-of-sight (LOS) path generation 
capability is formulated. In [Jun & He (2011)] an 
adaptive output feedback controller based on neural 
network feedback-feed forward compensator (NNFFC) 
which drives a surface ship at high speed to track a 
desired trajectory has beendesigned. A hybrid model 
was proposed in [Wu & Penga (2012)] to represent the 
ship's tracking dynamic behavior. A single-input 
single-output nonlinear time series model has 
beenbuilt to characterizethe responses between the 
ship's heading angle deviation and its rudder angle. 

Sliding mode control is based on the Lyapunov 
stability theorem and the LaSalles invariant principle, 
often referred toas variable structure.Sliding mode 
control (SMC) is a control method that switches 
between two distinctly different control laws 
depending on the state of the system. The salient 
advantages of sliding model control are: i) fast 
response and good transient performance; ii) its 
robustness against a large class of perturbations or 
model uncertainties and its ability to reject external 
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disturbances; and iii) the possibility of stabilizing 
some complex nonlinear systems which are difficult to 
stabilize by continuous state feedback laws. 

The SMC can be used as a supervisory controller along 
with the adaptive control. Adding the SMC to the 
control scheme will enable the overall controller to 
adapt over time as well as reject external disturbances 
and protect against parameter variations [Ebel (2011)], 
[Chan (2000)], [Lokukaluge & Perera (2012)]. This is of 
great importance in the case of ship control systems in 
which there are uncertainty of nonlinear model 
parameters in the presence of disturbances such as 
waves, wind and currents. 

The heading measurement is corrupted with colored 
noise due to wind, waves and ocean currents as well as 
sensor noise. However, only the slowly varying 
disturbances should be counteracted by the propulsion 
system, whereas the oscillatoric motion due to the waves 
should not enter the feedback loop. This is done by 
using so called wave filtering technique which separates 
the heading measurement into a low frequency and 
wave frequency estimate [Ricardo, Guedes (2009)]. For 
this purpose, the recursive filtersis utilized. 

In this paper, the problem of way-point tracking of 
container ship is considered. The disturbances such as 
wave, winds and currents are regarded. These 
disturbances are time-varying and control problem 
has been solved using recursive filters and adaptive 
stochastic sliding modemethod. 

The remainder of this paper is organized as follows. 
Section 2 deals with the modeling of ship and 
disturbances. Section 3 describes designing nonlinear 
controller, in section 4 the LOS guidance system is 
explained briefly. Section 5 discusses results and, 
finally, andthe conclusion is made in the last section. 

Equations of Motion and Disturbances 
Modelling 

Equations of Motion 

Ship dynamics is obtained by applying Newton's laws. 
The marine vehicle has 6 degrees-of-freedom (DOF) 
since six independent coordinates are necessary to 
determine the spatial position and orientation of a 
rigid body. The six different motion components are 
called: surge, sway, heave, roll, pitch, and yaw. 
Accordingly, the most generally used notation for 
these quantities ar e:x,y,z, cf>,dandip: "Fig. T'shows all 
six coordinate definitions and the most generally 
adopted reference frame. The position and orientation 


of the ship are described relative to the inertial 
reference frame 0 E - x E y E z E (Earth-fixed reference 
frame). The general ship equations of motion can be 
expressed in compact form as [Velagica et al. (2003)] 



FIG. 1 BODY- AND EARTH-FIXED REFERENCE FRAMES 


Mv + C(v) V = T (1) 

Where M is the inertia matrix, C(v ) is a matrix of 
Coriolis and centripetal terms, v — [u,v,w,p, q,r] T is 
the body-fixed (ship) linear and angular velocities 
vector and t = [ X , Y , Z, K, M, N] T is a generalized vector 
of external forces and moments. For the course- 
keeping and the track-keeping problems, only the 
horizontal-plane ship motion is used. 

Due to that, 6DOF model is simplified and reduced to 
3DOF model. Assuming that the coordinate of body- 
fixed frame origin is set in the center line of the ship 
(y G = 0), the mass distribution is homogeneous, the xz- 
plane is symmetrical, and that the influence of motions 
in the z-direction (heave), rotation about x-axis (roll) 
and y-axis (pitch) to the motion in horizontal plane can 
be neglected, the non-linear ship equations of motion 
are (Fossen, 1994): 


Surge: m(u — vr — x G r 2 ) = X 

(2) 

Sway: m{y + ur + x G f ) = Y 

(3) 

Yaw: I z r + mx G (v + ur) — N 

(4) 


Where m is mass of the ship, u and v are surge and 
sway velocities, respectively, r is the yaw rate, I z is 
moment of inertia about the z-axis, X and Y are forces 
in the x- and y-axis direction, respectively, N is a 
moment around the z-axis and r G — [x G ,y G ,z G ] is the 
center of the gravity. 

X = X(u, v,r, u, v,f, S) (5) 

Y = Y(u,v,r,u,v,f,S) (6) 

N = N(u, v, r, ii, v, r, S ) (7) 

Disturbances Modeling 

While moving in oceans, ships' motions are often 
influenced by environmental disturbances, therefore 
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in order to control ships effectively, it is necessary to 
model the environmental disturbances such aswinds, 
ocean currents and waves. Influence of each type of 
disturbances is derived separately and then the 
principle of superposition is applied to get the 
influences of environmental disturbances [Fossen 
(1994)]. 

1) Influences of Wind 

Forces and moments generated by the wind are 
given by the following system of equations: 


= (1/2 )p A c x ( e a )A T vf 

( 8 ) 

Y a = (1/2 )p A C y ( e a )A L vf 

( 9 ) 

N a = (1/2 )p A C n (6 a )LA L V 2 

(10) 


These forces and moments are added to right hand 
side of the "Equations (2) to (4)". 

Here, p A is air density; A T A L are transverse and 
longitudinal projected areas, respectively; 9 a , is 
wind relative direction; V A is relative wind speed 
and C x , C y , C n are forces and moment coefficients in 
X , Y , N directions, respectively. The relation 
between wind relative direction and forces and 
moment coefficients are shown in Fig. 2. 



FIG. 2 RELATION BETWEEN RELATIVE WIND DIRECTION AND 
FORCES AND MOMENT 

2) Ocean Currents 

The two-dimensional current model is used here. 
The Earth-fixed current components can be 
described by two parameters: the average current 
speed V c and direction of the currents y c . The body- 
fixed components can be computed from 

u c = V c cos(y c -xp) (11) 

v c = V c sin(y c - ip) (12) 

The average sea current velocity for computer 
simulations can be generated by using a first-order 
Gauss-Markov process, described by the following 


differential equation: 

Vc(t) + p Q V c (t) = w(t) (13) 

Where w(t) is a zero mean Gaussian white noise 
sequence and p Q > 0 is a constant. This process 
must be limited such that V min < V c (t') < V max in 
order to simulate realistic sea currents. 

3) Wind-Generated Wave 

The oscillatoric motion of the waves can be 
described by (Fossen, 1998): 

V'hOO = K w s/(s 2 + 2(a> 0 s + wg) X w H (s) (14) 

Where w H (s) is a zero-mean Gaussian white noise 
process, co 0 is a wave frequency (modal frequency), 
is a damping coefficient and K w is a gain constant. 
The gain constant is defined as k w = 2^'cl> 0 ct w ; where 
o w is a constant describing the wave intensity. 

The transfer function in "Equation (14)" is usually 
represented by the state space model 

■fw = h (15) 

= -2Sa> 0 i p H - WoG + K w w h 

Designing Nonlinear Controller 

Nonlinear Model of Ship 

The control system is designed for steering a ship on 
the course. In this system, the controlled parameter is 
the ship course, ip(t), while the controlling parameter 
is the rudder angle, 5(t). For control purpose, the 
Bech's model is used obtained from the second-order 
Nomoto model, in which the angular velocity ip(t) was 
replaced by a nonlinear manoeuvring characteristic 
H(xp(t)) r the coefficients of which can be determined 
from a spiral test that is used in many papers and 
described in (Fossen, 2000). The obtained model is 
given by the following equation (Amerongen, 1982): 

+ ^r) (^(O) 

= (^(t) +6(t)) (16) 

Where H( xp (t)) is approximated by the following 

function: 

H(ip(t))= oap(ty+ fixp(t) (17) 

Where a and fi are real constants. Commonly, a and fi 
are calculated based on a spiral test off-line, but in this 
article, because of the presence of time-varying 

disturbances, the spiral test is not helpful to obtain 
these parameters [Krestic et al. (1995)]. Identification 
of these parameters as well as the control procedure is 
implemented by adaptive stochastic sliding mode. To 
simplify control design, the Bech model is reduced to 
the Norrbin model given by: 
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TtKO + tf (<K0)=Ksa> (18) 

Where 

T=T ± +T 2 -T 3 (19) 

In order to carry out the system description by the 
nonlinear state equations, it is preferable to define the 
following state variables \x 1 — ip L (yaw angle), x 2 — 
xp — r L (yaw rate), and the output y =x 1 . The kinematic 
equations of ship dynamics are 


£l-^2 

(20) 

x 2 = a t x 2 + a 0 x 2 + du 

(21) 

u = b 

(22) 

The coefficient's values are defined as: 
% = —a/T 

(23) 

a 0 = —fd/T 

(24) 

d = K/T 

(25) 


Recursive Filters 

By combining the low frequency (LF) model 
representing the motion of the vessel and a high 
frequency (HF) lst-order wave induced motion, the 
heading angle can be expressed as: 

xp = xp L + xp H + v H (26) 


Wlute 

noise 



FIG. 3 LOW FREQUENCY (LF) AND HIGH FREQUENCY (HF) 
SUB MODELS 

Where v H represents zero-mean Gaussian measurement 
noise. The resulting model is shown in Fig. 3. In most 
ship applications, it is important that the contribution 
from the high frequency wave motion is suppressed. If 
not, wave disturbances can cause wear on the rudder, 
propeller and the thruster actuators. One of the 
methods for suppression of high frequency wave 
disturbances is state estimation. 

The ship wave system is described by the state 
x — VPl’ t L’^H’' 1 PhY ' input u = 6, and process noise 
w = [w l , w h ] t .It isassumed that w~lV(0, Q) and 
v~N(0, R). The nonlinear model is given by: 

x(n + 1) = <p[x(n),n\ + G[x(n),n]w(n) + B(n)u(n ) 
z(n) = h[x(n),n] + v(n ) (27) 


0 

0 


0 


T 

1 

0 

0 

0 

,B = 

K/T 

0 

,h = 

0 

0 

0 

K w . 


. 0 . 


. 1 . 


<p is the nonlinear function of combination ofthe 
"Equations (15), (20) and (21)". The state estimation in 
discrete nonlinear dynamic system has been appliedby 
recursive filters by the following equation (Schweppe, 
1973) 

x(N + 1\N + l)=x(N + 1|AQ + B(iV)u(A() + 

tf(lV + 1) X [z(N + 1) - h[x(N + 1\N),N + 1]} 
x(N + 1 1 A/) = <p[x(N\N), N] 

K(N + 1) = T(JV + 1| N + 1) 

X hM [x(N + 1\N),N + 1] R -1 (N + 1) 

T(JV + 1|7V + 1M H m [x(N + 1\N),N + 1] 

X IT 1 (N + 1 )[x(N + 1| N),N + 1] + + 1|A/)} _1 

Z(N + l|/V)=0 (1) [£(Al|Ai),Af]i:(A(|A()0 ( :h^ [x(N\N),N] 

+G [x (N | IV), N] Q (N)G ' [x (N \ A), N] 

T(0|0) = V 

£(0|0) = 0 (28) 
The control objectives are: 

(1) To force the output y=x 1 = xp of the system to 
asymptotically track the reference output y r (t) = xp r (t). 

(2) To keep the rudder angle in the acceptable range. 

Adaptive Stochastic Sliding Mode procedure 

The system is governed by stochastic differential 
equations in the Ito sense as 

dx 1 = x 2 dt 

dx 2 = [uiy| + a 0 x 2 ]dt + U{t)dt + dB(t) (29) 

Where 

£[dB(t)] = 0, 

E\dB(t,)dB(t 2 )\ = [ Wdt 0 ’ l=Jl = tz (30) 

We use the sliding function S = x 2 + A(x 1 — xp r ) (A > 
0). After obtaining the nominal sliding mode control, 
we focus on the following control 

f/(t) = — d 0 x 2 -a 1 x| — Ax 2 — kS (31) 

Where a 0 and a 1 are estimates of a 0 and a 0/ k > 0, and 
—kS is the 'switching' term. Instead of -G and x 2 ,xp(t) 
and f(t) which are obtained from state estimator 
should be replaced. 

Consider a Lyapunov function^ = E[V d ] where: 

+^(a° - a °) 2 - %) 2 (32) 

Clearly V > 0. To evaluate the time derivative of V d by 
following Ito lemma, we have 
1 1 

dV d = SdS + - (a 0 - a 0 )da 0 + - (a, - aAda-, 

Y V 

— S[(a 0 — cLq)x 2 + (a! — a 1 )x| — kS + 2 D]dt + 
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1 1 

S dB ^ t)t (a 0 CLo)dciQ E (u^ ir ^ ) (1. (i . i (33) 

Let 

da 0 = ySx 2 dt (34) 

= TjSx^dt (35) 

These define rules for updating parameter estimates 

a 0 and % . We have 

dV d = [-kS 2 + 2 D]dt + SdB(t) (36) 
Hence, we obtain 

V = dE[V d \/dt = —kE[S 2 ] + 2D (37) 


In the absence of stochastic excitation with D = 0, we 
have V < 0 . This shows that V is uniformly bounded. 

Los Guidance 

Systems for guidance consist of waypoints that are 
used to generate a trajectory (path) for the ship. A 
widely used method for path control is LOS guidance. 
In this methodology, a LOS vector from the ship to the 
next way point is computed. The desired heading 
angle as a set point for autopilot system can be 
calculated through: 

4> d (t) = tan-'Uy^k) -y(t))/(x d (k) -x(t))] (38) 

Care must be taken to select the proper quadrant for 
ip d . The next way point can be selected based on 
whether the vessel lies within a circle of acceptance 
with radius p 0 around the way-point (x d (k),y d (k) ) . 
Moreover, if the vessel location (*(t),y(t)) at the time 
t satisfies: 

[x d {k) - x(t)} 2 + \y d (k) - y(t)] 2 < Po 2 (39) 

The next way point (x d (fc + l),yd(k + 1)) should be selected. A 
guideline could be to choose p 0 equal to two shiplengths, that is 
p 0 = 2L [Fossen (1998)]. 

Simulation Results and Analysis 

Simulation results are considered for a container ship 
that the values of the parameters used during 
simulation are chosen from example 1.3 injFossen 
(1994)] and shown in Figs 4-9. For the numerical 
simulations, the fourth-order Runge-Kutta method 
with time-step h = 10 -4 is used. Total heading angle is 
shown in Fig. 4 that is obtained by stochastic sliding 
mode control without using recursive filters. The 
comparisons of LF heading angle with two control 
methods that are obtained from state estimation are 
shown in Fig. 5. Fig. 6 and Fig. 7 show the rudder 
angle of ship in the presence of time-varying 
disturbances by spiral test + PID method and adaptive 
stochastic sliding mode respectively, also the recursive 
filter is applied to both of them. Fig. 8 gives the way- 
point tracking course for the adaptive stochastic 
sliding mode autopilot with recursive filter in the 


presence of disturbances and without disturbances. 
These two courses are close to each other. Fig. 9 shows 
the updating parameters in adaptive stochastic sliding 
mode control. 



FIG. 4 COMBINING THE LOW FREQUENCY (LF) MODEL 
REPRESENTING THE MOTION OF THE VESSEL AND A HIGH 
FREQUENCY (HF) WAVE INDUCED MOTION. 



FIG. 5 COMPARISON OF HEADING ANGLE (LF) : SPIRAL TEST 
FOR OBTAINING SHIP MODEL PARAMETERS + PID 
CONTROLLER (DASHED LINE) AND ADAPTIVE STOCHASTIC 
SLIDING MODE(SOLID LINE) 


0.2 



- 0.2 1 1 1 1 1 1 

0 500 1000 1500 2000 2500 3000 3500 

Time(s) 


FIG. 6 SHIP RUDDER ANGLE IN THE PRESENCE OF 
DISTURBANCES BY SPIRAL TEST +PID 
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FIG. 7 SHIP RUDDER ANGLE IN THE PRESENCE OF 
DISTURBANCES BY ADAPTIVE STOCHASTIC SLIDING MODE. 



FIG. 8 COMPARISON OF SIMULATION RESULTS OF WAY- 
POINT TRACKING : IN PRESENCE OF DISTURBANCES 
(DASHED LINE), WITHOUT DISTURBANCES (SOLID LINE). 



FIG. 9 UPDATING THE PARAMETERS a 1; a 0 


Conclusions 

The paper has considered way-point tracking system in 
the presence of time-varying disturbances with 
adaptive stochastic sliding mode approach. Also for 


suppression of high frequency of rudder angle, the 
recursive filter is used. Even though PID control is 
widely used throughout industrial application, it 
handles parameter variations and external disturbances 
poorly, whereas the adaptive stochastic sliding mode 
control is known for its robustness against these 
uncertainties. The maximum distance of way-point 
tracking course for the adaptive stochastic sliding 
mode autopilot in the presence of disturbances and 
without disturbances is 300 m that is acceptable for 
this container ship with length 175 m.To compare the 
adaptive stochastic sliding mode control with PID, the 
recursive filter is applied to both of them. The results 
show the better performance in terms of transient and 
saturation in the control input for the adaptive 
stochastic sliding mode in comparison with PID (that 
the model parameters are obtained off-line by spiral 
test) because the model parameters are updated 
duringprocedure. Maximum rudder angle for this 
container ship is 10 (deg) that in adaptive stochastic 
sliding mode control, the rudder angle is saturated in 
the shorter time compared to PID control. So this 
algorithm iseffectively implementedin the way-point 
tracking control. 
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